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Ab-initio description of the magnetic shape anisotropy due to the Breit interaction 
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A quantum-mechanical description of the magnetic shape anisotropy, that is usually ascribed to 
the classical magnetic dipole-dipole interaction, has been developed. This is achieved by including 
the Breit-interaction, that can be seen as an electronic current-current interaction in addition to 
the conventional Coulomb interaction, within fully relativistic band structure calculations. The 
major sources of the magnetic anisotropy, spin-orbit coupling and the Breit-interaction, are treated 
coherently this way. This seems to be especially important for layered systems for which often both 
sources contribute with opposite sign to the magnetic anisotropy energy. Applications to layered 
transition metal systems are presented to demonstrate the implications of this new approach in 
treating the magnetic shape anisotropy. 
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Magnetic anisotropy is among the most important 
properties of magnetic materials in particular concern- 
ing their application in devices. When discussing the 
magnetic anisotropy energy of a material, denoting the 
difference in energy for two orientations of the magneti- 
sation, an incoherent approach is used so far [l], [J. On 
the one hand side, the dependency of the electronic struc- 
ture and the associated total energy on the orientation 
of the magnetisation, that is induced by spin-orbit cou- 
pling (SOC), is accounted for by corresponding relativis- 
tic band structure calculations. On the other hand, the 
additional shape anisotropy is ascribed to the anisotropy 
of the magnetic dipole-dipole coupling, that is treated in 
a classical way. This hybrid approach is used in particu- 
lar when dealing with layered transition metal systems. 
The pioneering and successful theoretical work on mag- 
netic surface films by Gay and Richter [3| was followed 
later on by many more investigations that benefited from 
the extension of standard band structure schemes to ac- 
count simultaneously for the presence of spin-orbit cou- 
pling and spin-polarisation, i.e. magnetisation, in a nu- 
merically reliable way. Especially interesting in this con- 
text are investigations on systems showing a competi- 
tion of the SOC and shape induced contributions to the 
magnetic anisotropy energy. This situation is frequently 
encountered for magnetic multi-layer or surface layer sys- 
tems for which a flip of the magnetic easy axis from out- 
of-plane to in-plane may be observed when the thickness 
of the magnetic layer is increased starting from a single 
mono-layer. In fact corresponding experimental findings 
could be reproduced by calculations based on the above 
mentioned hybrid scheme in the case of magnetic multi- 
layers [H as well as surface layer systems 

In spite of the successful applications of this hybrid 
scheme in dealing with the magnetic anisotropy, one has 
to keep in mind that there is no real justification for its 
use. Furthermore, as there have been no coherent quan- 
tum mechanical investigations performed so far, there is 
no experience on the range of its applicability. It seems 



that the first steps towards a coherent quantum mechan- 
ical description of the magnetic anisotropy were made by 
Jansen @, 0] , who pointed out that the shape anisotropy 
is ultimately caused by the Breit interaction - a relativis- 
tic correction to the Coulomb interaction between moving 
electrons [1, 0] . While Jansen performed model calcula- 
tions to investigate the magnetic anisotropy energy con- 
tribution caused by the Breit interaction, first numerical 
investigations were done by Stiles et al. [loj . However, 
these authors restricted to the spin-other-orbit part of the 
Breit interaction that implies a current-current interac- 
tion. Performing a Gordon decomposition of the current 
density into a spin and orbital part 11( one gets three 
additional terms. Moreover, Stiles et al. investigated the 
pure bulk ferromagnetic metals bcc-Fe, hep-Co and fcc- 
Ni for which the shape anisotropy is rather small. As a 
consequence, the results of this first numerical work are 
not very conclusive. 

In this contribution we present a fully relativistic de- 
scription of the magnetic anisotropy that accounts for 
the shape anisotropy by incorporating the full Breit in- 
teraction within the Dirac equation for magnetic solids. 
The scheme has been implemented by using a fully rel- 
ativistic multiple-scattering or Korringa-Kohn-Rostoker 
(KKR) formalism. To demonstrate the power of this new 
approach applications to the layered systems Fe„Pd ra and 
Fe n /Au(001) are presented. The results for the magnetic 
anisotropy energy are compared to results obtained using 
the conventional hybrid approach. 

The Breit-interaction ||, Q is a correction to the 
Coulomb interaction between moving electrons that may 
be split into its magnetic and retardation part: 
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Within this relativistic formulation, the standard Dirac 
matrices cvi are connected to the current density operator 
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ji via j = — \e\ca [ll| and tn denotes the distance be- 
tween two electrons. The Gordon decomposition of the 
electronic current mentioned above is not used in the fol- 
lowing. This implies that all coupling mechanisms con- 
nected with the Breit interaction are finally accounted 
for. 

In the following we focus here on the magnetic part as 
the retardation term docs not contribute on the Hartree 



level [12j. This allows to represent the interaction of an 
electron with all the other electrons by the vector poten- 
tial: 



Mr ) = l ~ I rfV^l 
c \r — r 



(2) 



Including the corresponding interaction term within elec- 
tronic structure calculations done in the framework of rel- 
ativistic spin density functional theory (SDFT) [l3[ leads 
to the following single particle Dirac Hamiltonian: 

H D = ~icd-V + /3mc 2 + V c ff(r) + (3Z-B(r) + ed-A(r) , 

(3) 

with the corresponding total energy E given by: 

E = T s + E cxt + £g + El + E xc (4) 

Here the first three and the last terms have their usual 
meaning, i.e. they stand for the kinetic energy (T s ), 
the electron- nucleus (E CKt ) and electron-electron 
Coulomb interaction and the exchange-correlation (E xc ) 
contributions. The fourth term: 
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going beyond standard SDFT, is a Hartree-like term due 
to the Breit or current-current interaction. As mentioned 
above the retardation term in Eq. ([!} does not contribute 
to El pj| . This could indeed be veryfied by our numer- 
ical results. 

As usually done for total energy calculations, the ki- 
netic energy term in Eq. (jj) can be eliminated by making 
use of the Hamiltonian in Eq. ©. 

When calculating the vector potential A(f ) it is advan- 
tageous to decompose the integration regime in Eq. @ 
into the central atomic cell i and its surrounding giving 
rise to the on- and off-site contribution, respectively, for 
the vector potential within cell i: 



A\{r) = AT{r) + Af(f). 



(6) 



The on-site contribution A° n can be determined directly 
from the currents within atomic cell i, on the basis of 
Eq. ((T|). The off-site contribution A° s from all other 
sites can be obtained by applying the common far field 
approximation 



Af(r) = 



(7) 



where the total magnetic moment Mj represents the cur- 
rent distribution in atomic cell j. The lattice summation 
in Eq. ([7]) is dealt with by an Ewald summation tech- 
nique in the case of two- or three-dimensional periodic 
systems. 

Using a band structure method based on a decompo- 
sition of the system into atomic cells, as for the KKR- 
method used here, Eq. ([3]) has to be solved in a first step 
for isolated atomic cells (single-site problem). For this 
purpose it is helpful to expand the vector potential by 
means of spherical harmonics. For the Breit-interaction 
term, the last term in Eq. one obtains: 



n 



BI 



e ^2 a r, 



E 

l A m A 



A-m a 



(r)Y™*(r) 



(8) 



In the implementation presented here the atomic 
sphere approximation (ASA) together with a restriction 
to collinear magnetism has been applied, i.e. within an 
atomic cell one has V e s(r ) — V(r) and B(r ) = B(r)eM, 
with the orientation of the magnetic moments ejj that 
fixes the local z-axis. In line with these geometrical sim- 
plifications A(f ) is restricted to have rotational symme- 
try around eu with A(r) pointing everywhere in tangen- 
tial direction, i.e. A(r) = A(r,8)e ( f > implying A^_ x {r) = 
— Aj^ +1 (r) and all other terms being zero. (Further tech- 
nical details will be given elsewhere [14j). 

With the single-site Dirac equation being solved, the 
electronic structure of the investigated system can be 
calculated by means of the standard relativistic multi- 
ple scattering or KKR-method [HI]. Calculating the to- 
tal energy E(h) on the basis of Eq. Q for two different 
orientations of the magnetisation n and h' , respectively, 
gives the corresponding magnetic anisotropy energy AE 
as the difference E(n) — E(n'). 

As an application of the approach sketched above we 
consider the geometrically simple case of a free-standing 
bec Fe monolayer with a lattice constant of bulk Fe 
(a = 2.87 A). The magnetic moment M was taken to 
point out-of-plane, i.e. along the z-axis of the system. 
The top left panel of Fig. Q] shows the resulting radial 
vector functions AjJ^ A (r) for tua — 1 with the on- and 

off-site parts to A(r) indicated by black and red lines, 
respectively, 

The top right panel of Fig. [T] shows the resulting radial 
dependence of the ^-component A™ along the z-direction 
(black line) together with A°^ obtained via Eq. ([7]) (red 
line). Applying the far field approximation of Eq. 
also to the on-site term results in the blue curve which 
is divergent at the origin. In the outermost region of an 
atomic sphere, however, the far field approximation to 
A(f) is already very good justifying the use of Eq. (0 
also for nearest neighbouring atoms. 

In order to verify the correctness of A™ when calcu- 
lated directly via Wbi (see Eq. ([T])) we also determined 
the current density j(f) within an Fe sphere and then 
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FIG. 1: Vector potential and current density for Fe in a free- 
standing monolayer. Top left: On- (black) and off-site (red) 
contributions to the vector potential functions AJ^^ A (r) for 
I a = 1,3,5 and tua — 1 in atomic units. Top right: On- 
(black) and off-site (red) component along the x-direction 
[6 = ^,<f} = 0). A on is compared to the far field approxima- 
tion A s * (blue curve). Bottom left: radial current density 
functions for i — 1,3,5 and m = 1. Bottom right: 

comparison of A$ along x calculated via Hbi (full line) and 
from the current density (dashed line). 



computed A™ via Eq. ((2]). The bottom left panel of 
Fig. [1] shows the radial electronic current density distri- 
bution functions jj^if) (defined in analogy to Aj™£ 
in Eq. ((5])) for an Fe atom and the bottom right panel 
of Fig. [T] presents a comparison of A™ resulting from 
jim( r ) (dashed line) and obtained from Hbi- As one 
can see both approaches give nearly identical results. 

A further impression of the spatial variation of A(f) 
and j(r) is given in Fig. [5] showing the vector fields in 
the xy-plane and their color-coded amplitude in the xz- 
plane. The figure reflects the rotational symmetry of 
A(r) imposed by the use of the ASA (atomic sphere 
approximation) and the alignment of the magnetisation 
along the z-direction. 

For the anisotropy energy AE = E(x) — E(z) of the 
free-standing bcc Fe monolayer we obtain —0.063 meV 
with the sign indicating that the magnetisation favours 
an in-plane orientation. AE can be decomposed into 
AEsoc = 0.096 meV preferring an out-of-plane mag- 
netic easy axis and a dominating AE-qi — —0.159 meV 
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FIG. 2: Top: Vector potential A(f) of Fe within the xy-plane 
with z = (left) and within the :rz-plane with y = (right) for 
the magnetisation along the 2-axis. Bottom: Corresponding 
current density j(r) within the a;j/-plane (left) and within the 
xz-plane (right). 



causing the magnetisation to lie in-plane. Here, AEbi = 
AE — AEsoc has been obtained by performing calcula- 
tions with (SOC+BI) and without (SOC only) Hbi giving 
AE and Ai?soc, respectively. For the classical dipolar 
shape anisotropy AE^d we get —0.169 meV which agrees 
astonishingly well with the quantum mechanical A£bi- 

A second application dealt with the multilayer systems 
Fe„Pd„. For these a fee structure with (OOl)-oriented 
atomic layers has been assumed, i.e. with the z-axis per- 
pendicular to the Fe- and Pd-layers, respectively. For 
the special case n = 1 this corresponds to the CuAu- 
structure. 

The impact of the Breit-interaction on the electronic 
structure leads to a competition of the various contribu- 
tions to the magnetic anisotropy energy AE of Fe n Pd n 
when the parameter n is varied. Using the conventional 
approach that accounts only for AEsoc within the elec- 
tronic structure calculations leads to a strong contribu- 
tion that favors an out-of-plane orientation of the mag- 
netisation for n = 1 — 6 (see Fig. [3]). The additional 
classical shape anisotropy contribution AE^d gives rise 
to a contribution that favours an in-plane orientation 
of the magnetisation. As one may expect, AEdd of 
Fe n Pd n is primarily determined by the magnetic mo- 
ments within the Fe layers that amounts to be between 
2.74 and 2.95 fiB for all values of n considered here. The 
induced moments on the Pd layers that are in the range 
0.005 to 0.330 hb are only of minor importance. As a re- 
sult of this AEdd increases nearly linearly with n for the 
range of n considered here. For n — 5 and presumably 
also for n > 6 the magnetic dipole-dipole term exceeds 
the SOC-induced term leading to a flip of the magnetic 
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FIG. 3: Magnetic anisotropy energies for repeating Fe„Pd n 
multilayers as function of n: total magnetic anisotropy energy 
AE (squares) and its decomposition into magneto-crystalline 
part A£soc (circles) and Breit part AEbi (diamonds) which 
is compared to its classical approximation (crosses). 
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tern Fe„/Au(001) that exhibits a flip of the easy axis from 
out-of-plane to in-plane for around three monolayers of 
Fe (see also Ref. [5|). Again the contribution AEbi to 
the total magnetic anisotropy energy AE is found to be 
very close to the classical result AE^d- Only for rather 
large thicknesses a noteworthy deviation of the two can 
be seen in Fig. [4] 

In summary the Breit interaction has been incorpo- 
rated within fully relativistic band structure calculations 
for magnetic layered systems. This development gives 
access to an ab-initio calculation of the magnetic shape 
anisotropy energy using a coherent approach that ac- 
counts simultaneously for spin-orbit coupling and the 
Breit interaction. First applications of this new approach 
to the systems Fe„Pd„ and Fe„/Au(001) were pre- 
sented. Taking the difference of the calculated magnetic 
anisotropy energy obtained via the combined SOC+BI - 
and the SOC-only approaches, the contribution due to 
the Breit interaction could be separated. For all systems 
investigated so far it was found that the resulting Breit 
contribution is very close to the classical result calculated 
on the basis of the magnetic dipole-dipole interaction. 
This result can be explained to some extend by the fact 
that it is the magnetic part of the Breit interaction that 
gives rise to the shape anisotropy while the retardation 
term does not contribute on the Hartree-level. 
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gemeinschaft) project EB 154/14-2 is gratefully acknowl- 
edged. 



FIG. 4: As for Fig. [3] but for the surface layer sytem 
Fe n /Au(001). 



easy axis from out-of-plane to in-plane. 

Calculating the anisotropy energy AE on the basis of 
the coherent SOC+BI scheme, AE follows qualitatively 
the variation of the SOC-induced magnetic anisotropy 
energy with n and implies also a flip of the magnetic 
easy axis (see Fig. [3]). To allow for a direct comparison 
of the two approaches the difference AEbi of SOC+BI 
and the SOC-only scheme is shown as well in Fig. [H As 
can be seen, AEbi is very close to the classical AE^d- 
This result obviously justifies the use of the conventional 
classical approach for the shape anisotropy used so far. 
In particular the conventional scheme seems to reproduce 
the quantum mechanical result not only in a qualitative 
but in general also quantitatively in a satisfying way. Ob- 
viously, only for rather short interatomic distances one 
has to be aware of possibly pronounced deviations be- 
tween the classical and quantum mechanical approaches. 
These conclusions are confirmed by results obtained for 
the closely related multilayer systems Fe„Pt„, Co„Pd n 
and Co„Pt„ as well as first applications to surface layer 
systems. Fig. U shows corresponding results for the sys- 
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